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1. Introduction 

Quantum graphs attracted a lot of interest recently. There are several rea- 
sons for that. On one hand these models are useful as descriptions of vari- 
ous structures prepared from semiconductor wires, carbon nanotubes, and 
other substances. On the other hand they provide a tool to study proper- 
ties of quantum dynamics in situations when the system has a nontrivial 
geometrical or topological structure. 

Quantum graph models contain typically free parameters related to cou- 
pling of the wave functions at the graph vertices, and to get full grasp of the 
theory one has to understand their physical meaning. A natural approach 
to this question is to investigate "fat graphs" , that is, systems of thin tubes 
built over the skeleton of a given graph, and to analyze its limit as the tube 
thickness tends to zero. 

While simple at a glance, the problem is in fact rather difficult and its 
understanding is being reached through a long series of works. The aim of 
the present paper is to review some recent achievements in this area. We 
present this survey in a non-technical way referring for detailed proofs and 
a wider background to the literature, in particular, to our recent papers 
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[7,19]. Having said that it is important to stress that we will formulate the 
problem and the results in a fully rigorous way. 



2. Quantum graphs 
2.1. A bit of history 

The quantum graph concept was born in early days of quantum mechanics 
being hrst suggested in the 1930's by Linus Pauling as a model of aromatic 
hydrocarbons, and worked out later by Ruedenberg and Scherr [37]. Then, 
as it sometimes happen in the history of science, it was happily forgotten. 

In a sense it might be surprising because the idea of a quantum particle 
living on a graph is theoretically attractive, however, it was not enough 
and for three decades quantum graph models enjoyed the status of an ob- 
scure textbook example. This changed in the eighties, when the diminishing 
size of structures produced in solid-state-physics laboratories reached the 
state when the electron transport in them became dominantly ballistic and 
quantum graphs suddenly reemerged as a useful model. 

The list of physical system to which these methods can be applied kept 
expanding. At the beginning it included microstructures fabricated from 
semiconductor or metallic materials, later carbon nanotubes were added. 
It is worth mentioning, however, that the same technique can be used also 
to investigation of electromagnetic phenomena in large network-type struc- 
tures [27], at least as long stationary situation is considered. 

Quantum dynamics of a particle confined to a graph can mean var- 
ious things, of course. Typically one considers a nonrelativistic situation 
described by a Schrddinger operator supported by the graph. Often the 
motion is free but in other situations one adds potentials corresponding to 
external electric or magnetic fields, spin degrees of freedom, etc. Graphs 
can support also Dirac operators. Such a model, too, was for a long time 
regarded as a theoretician toy and attracted a limited attention only [3,4]. 
The situation changed dramatically two or three years ago with the discov- 
ery of graphene in which electron behave effectively as relativistic particles 
which triggered a wave of papers of the subject. 

The literature on quantum graphs is nowadays immense; we are not 
going to try to give a bibliographical review and refer instead to the pro- 
ceedings volume of a recent Isaac Newton Institute programme [15] where 
one can find an extensive guide to further reading. 
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2.2. Vertex coupling 

For simplicity let us consider a graph having a single vertex in form of a 
star, i.e. n halfiines with the endpoint connected. The state Hilbert space 
H of such a system is = 0™ =1 L 2 (R + ) and the particle Hamiltonian acts 
on W as ipj n- —ip"] the values of physical constants are irrelevant for our 
discussion and we put conventionally h = 2m = 1 . 

The Hamiltonian domain consists of W 1 ' 2 functions; in order to make it 
self-adjoint we need to impose suitable boundary conditions at the vertex. 
Since we deal with a second-order operator, the latter involve boundary 
value ^(0) := {^>j(0)} and ^'(O) := {tpj(0)}; conventionally they are writ- 
ten in the form 

^*(0) + B*'(0) =0 (1) 

proposed by Kostrykin and Schrader [28], where the n x n matrices A, B 
give rise to a self-adjoint operator iff they satisfy the conditions 

• rank (^1|-B) = n 

• AB* is self-adjoint 

The obvious drawback of (1) is that the pair A, B is not unique. The com- 
mon way to remove the non-uniqueness [24,26,29] is to choose 

A = U-I, B = i(U + I), (2) 

where U is an n x n unitary matrix; there are also other unique forms 
more suitable for some purposes [7,8,30]; one of them we will need in Sec. 6 
below. It is obvious from (2) that the coupling of n edges is characterized 
in general by n 2 real parameters. 

There is a simple way to derive the boundary conditions with which 
can be traced to [23] where it was used for n = 2. Self-adjointness requires 
vanishing of the boundary form, 

n 

£(^-$^)(o) = o, 

which occurs iff the norms H^O) ± ^^'(O^lc™ with a fixed t ^ coincide, 
since the difference of the squared norms is just the Ihs of the displayed rela- 
tion. Consequently, the vectors must be related by an n x n unitary matrix, 
which yields immediately (U - /) V E'(0) + i£(U + 7)*'(0) = 0. It may seem 
that we have an extra parameter here, however, matrices corresponding to 
two different values of i are related by 

_ {i + (')U + t-i' 

~ {£ - l')U + 1 + £' ' 
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so it just fixes the length scale of the problem and we can put 1=1 
without loss of generality. Note also that the parameter matrix is closely 
related to the scattering at the vertex, specifically, it coincides with the 
on-shell scattering matrix at the momentum k = 1. 

2.3. Examples of vertex coupling 

Denote by J the n x n matrix whose all entries are equal to one; then 
the unitary matrix U — n + ia J — I corresponds to the standard S coupling 
characterized by the conditions 

n 

^•(0)=V*(0)=:V(0), j,k = l,...,n, ^V;(0)=<#(0) (3) 

j=i 

of "coupling strength" a G R; we include also the case a — oo, or U = —I, 
when the edges are decoupled with Dirichlet conditions at the endpoints. 
Another particular case of interest is a — corresponding to the "free 
motion". It would be natural to call then (3) free boundary conditions, 
however, they are mostly called Kirchhoff in the literature a . Note that the 
f5-couplings are the only ones with wave functions continuous at the vertex. 

The second example to mention is the S' s coupling, a counterpart to the 
above one with the roles of functions and derivatives interchanged. The 
corresponding unitary matrix is U = I — n _^ >7 giving 

n 

^.(0)=^(0)=:V'(0), j,k = l,...,n, 5>i(°)=^'(°) (4) 

3=1 

with P G R; for j3 ~ oo we get decoupled edges with Neumann conditions. 

3. Vertex understanding through approximations 
3.1. Statement of the problem 

The first question to pose is why we should be interested in quantum graph 
vertex couplings. There are several reasons for that: 

• One is mathematical. Different couplings define different Hamilto- 
nians which have different spectral properties. Sometimes they can 
be quite involved; as an example let us number theoretic properties 
of rectangular lattice-graph spectra [13]. 



a The name is generally accepted but unfortunate because in electricity it is associated 
with current conservation at the junction, and in the quantum case any self-adjoint 
coupling preserves probability current. 
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• On a more practical side, the conductivity of nanostructures is 
controlled typically by application of external fields. Understanding 
of vertex coupling would give us an alternative mean to this goal. 

• As a specific example, the authors of Ref. 10 used the generalized 
point interaction on line as a model of a qubit; in a similar way 
star graphs with n > 2 edges can similarly model qudits. 

At a glance the vertex parameters can be interpreted easily. One should 
replace the graph in question by a family of "fat graphs" , i.e. a tube net- 
work built around the graph skeleton, with appropriate Laplacian as the 
Hamiltonian. Such a system has no free parameters, so it would be enough 
to inspect the squeezing limit with the tube diameter tends to zero and to 
see which graph Hamiltonian we obtain. Unfortunately, as it is often the 
case with simple answers, the problem is in reality rather complicated: 

• The answer depends substantially on the type of the Laplacian sup- 
ported by the tube network. The Neumann case is easier and after 
an effort more than a decade long an understanding was reached 
[17,18,22,31,34,36,38]. The drawback was that the limit gave the 
free (Kirchhoff ) boundary conditions only. 

• the Dirichlet case is more difficult and only recently some substan- 
tial results were obtained [1,5,12,25,32,33], nevertheless, a lot of 
work remains to be done 

Before proceeding to our main topic, let us review briefly the existing results 
we have mentioned above. 

3.2. Briefly on Dirichlet networks 

The distinctive feature of the Dirichlet case is the energy blow-up associ- 
ated with the fact that the transverse part of the Dirichlet Laplacian has 
lowest eigenvalue proportional to d~ 2 where d is the tube diameter. To get 
a meaningful result we have thus to use an energy renormalization which 
can be done in different ways. Molchanov and Vainberg [32] chose the en- 




Fig. 1. Energy renormalization 
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ergy A to be subtracted between the first and second transverse eigenvalue, 
cf. Fig. 1, and obtained a nontrivial limit determined by scattering prop- 
erties of the corresponding "fat star" . A drawback of this approach is that 
leads to energy spectrum unbounded from below which is a feature one tries 
to avoid in meaningful nonrelativistic models. 

Most authors choose therefore the transverse threshold Ai as the energy 
to subtract. In such a case the limit is generically trivial giving disconnected 
edges with Dirichlet endpoints [32,33]. However, the limit can be nontrivial 
provided the tube system we start with has a threshold resonance [1,5,25]; 
a similar, closely related effect using finite star graphs was proposed in [12]. 

3.3. A survey on Neumann network results 

Consider first for simplicity a finite connected graph Mq with vertices Vk, 
k e K and edges ej ~ Ij :— [0, ij], j € J; the corresponding Hilbert space is 
thus L 2 (M ) := 0, ,/.-:/,<• The form u ^ \\u'\\l Io ■= J2 jeJ with 
u £ W 1,1 (Mq) is associated with the operator which acts as —Am u = —u" 
and satisfies the free boundary conditions. 




Fig. 2. Graph Mo and fat graph M e 



On the other hand, consider a Ricmannian manifold X of dimension 
d > 2 and the corresponding space L 2 (X) w.r.t. volume dX equal to 
(det g) l / 2 dx in a fixed chart. For u £ C^ mp (X) we set 

qx (u) := \\dufx = I |du| 2 dX, |du| 2 = S^gVdiudju; (5) 

the closure of this form is associated with the self-adjoint Neumann Lapla- 
cian Ax on the X. Let us stress that within this framework we can treat 
both "solid" tubes with the boundary at which Neumann condition is im- 
posed, as well as "sleeve-type" manifolds without a boundary when the 
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particle is supposed to live on the surface - cf. Fig. 2. This is made possible 
by the similarity of the transverse ground-state eigenfunction in both cases. 

The "fat graphs" M £ associated with the graph M are all constructed 
from X by taking a suitable e-dependent family of metrics. This the ap- 
proach was used in [17]; in contrast to earlier work such as [31] one also 
need not assume that the network is embedded in a Euclidean space since 
only intrinsic geometrical properties are involved. 

The analysis requires dissection of M e into a union of edge and vertex 
components, U £ j and V Sj k, respectively, with appropriate scaling properties, 

• for edge regions we assume that U £ j is diffcomorphic to Ij x F 
where F is a compact and connected manifold (with or without a 
boundary) of dimension m := d — 1 

• for vertex regions we assume that the manifold V e ,k is diffcomorphic 
to an e- independent manifold Vk 

In this setting one can prove the following result [17]: 

Theorem 3.1. Under the stated assumptions we have eigenvalue conver- 
gence, \k{M e ) — > Afe(M ) , k = 1, 2, . . . , as e — V 0. 

The shrinking limit thus leads to free boundary conditions only, but also 
in other respects the stated result is not particularly strong, for instance, 
in that it concerns the eigenvalue convergence in finite graphs only. One 
can do better: in Ref. 34 Olaf Post proved a norm-resolvent convergence 
Am e Am as e — > 0+ on generally infinite graphs under natural uni- 
formity conditions analogous to those of used in Theorem 4.3, namely (i) 
existence of nontrivial bounds on vertex degrees and volumes, edge lengths, 
and the second Neumann eigenvalues at vertices, (ii) appropriate scaling 
(analogous to the described above) of the metrics at the edges and ver- 
tices. The involved operators act on different Hilbert spaces, of course, and 
the stated limiting relation makes sense with a suitable identification map 
which we will describe below. 

Other extensions are possible. For graphs with semi-infinite "outer" 
edges, e.g., the problem typically exhibits series of resonances, and one 
may ask what happens with them if the graph is replaced by a family of 
"fat" graphs. Using exterior complex scaling in the "longitudinal" variable 
one can prove a convergence result for resonances as e — > [18]; the same 
is true for embedded eigenvalues of the graph Laplacian which may remain 
embedded or become resonances for e > 0. 

Hence we have a number of convergence results is available for squeez- 
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ing limit of Neumann-type thin tube networks, however, the limiting oper- 
ator corresponds always to free boundary conditions only. The question is 
whether one can do better. 

4. Beyond the free coupling 

4.1. A graph inspiration 

It is obvious that one has to add a new feature to the approximating family 
to get more general results. Let us look how one can approximate d coupling 
on graphs using families of scaled potentials. For simplicity we will consider 
again the n-edge star graph as in Sec. 2.2, however, we replace the Laplacian 
at the edges by a Schrddinger operator, tpj i-> —ipj+Vjipj. In order make the 
problem well-defined we have to impose requirement on the potential; we 
suppose that Vj £ £ 1 1 oc (R+) , j = 1, . . . , n. If the boundary conditions at the 
vertex are (3) with the parameter a £ R we get a self-adjoint operator which 
we denote as H a {V). Let now the potential contain a scaled component, 

W*j-=\Wi(^) > 3 = l,...,n, (6) 
then we have the following result [14]: 

Theorem 4.1. Suppose that the potentials Vj £ i 1 1 oc (R+) are below 
bounded and Wj £ L 1 (R + ) for j = 1, . . . , n . Then 

H (V + W e )^H a (V) 

as e — >• 0+ in the norm resolvent sense, with the coupling parameter defined 
asa:= £"=i Wj (x) dx. 

Our aim is to "lift" this result to tube networks. 

4.2. Single vertex networks 

Consider first a star graph again. Let G have one vertex v and degv 
adjacent edges of lengths £ e £ (0, oo]. The corresponding Hilbert space 
is L 2 (G) := ® e££; L 2 (/ e ), the decoupled Sobolev space of order k is 
W'maxCG) : = ®eeE W 2 ' k (I e ) together with its natural norm. 

Let p = {p e } have components p e > for e £ E; we introduce it because 
we want to consider squeezing limits also in the situation when the tubes 
have different cross sections. The Sobolev space associated with weight p is 

W 2 p \G) :={/eW^(G): / £ Cp} , 
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where / := {/ e (0)}, in particular, if all the components are equal, p = 
(1, . . . , 1), we arrive at the continuous Sobolev space W 2 ' 1 (G) := W^' 1 {G). 

Next we have to introduce operators on the graph. We start with the 
(weighted) free Hamiltonian Aq defined via the quadratic form = J>g, 

Kf)--=\\f\\ 2 G = J2\\ e '\\l and doma^W^G) 

e 

for a fixed p (we drop the index p); the form is a closed as related to the 
Sobolev norm ||/|| wn,i(G) = II/'IIg+II/IIg- The Hamiltonian with S- coupling 
of strength q is defined via the quadratic form rj = f)(G, g ) given by 

W)--=\\f\\G + <l(v)\m\ 2 and domi,:=W^(G) 

Using standard Sobolev arguments one can show [19] that the 5-coupling 
is a "small" perturbation of Aq by estimating the difference f)(/) — f(/). 

The manifold model of the "fat" graph is constructed as in the previous 
section. Given e G (0, £o] we associate a d-dimensional manifold X E to the 
graph G in the following way: to the edge e G E and the vertex v we ascribe 
the Riemannian manifolds 

X £ , e := I e x sY e and X £iV := eX v , 

respectively, where eY e is a manifold Y e equipped with metric h e ^ e :— e 2 h e 
and eX £yV carries the metric g EtV = e 2 g v . As before, we use the e- 
independent coordinates (s,y) G X e = I e x Y e and x G X v , so the radius- 
type parameter e only enters via the Riemannian metric. Let us stress this 
includes the case of the e- neighborhood of an embedded graph G C M d , but 
only up to a longitudinal error of order of e; this problem can be dealt with 
again using an e-dependence of the metric in the longitudinal direction. 
The Hilbert space of the manifold model L 2 (X e ) can be decomposed as 

L 2 {X £ ) = (L 2 (I e ) ® L 2 (eY e )) (B L 2 (eX v ) 

e 

with the norm given accordingly by 

\W\\x e - E ^ I M 2d ^ds + e d f \u\ 2 dx v , 

eeE 

where dx e — dy e ds and dx v denote the Riemannian volume measures asso- 
ciated to the (unsealed) manifolds X e = I e x Y e and X v , respectively. We 
also introduce the Sobolev space W 1,1 {X e ) of order one defined conven- 
tionally as the completion of the space of smooth functions with compact 
support under the norm IM!^ w x , = ||du||^- + . 
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Next we pass to operators on the manifold. The Laplacian Ax e on X £ 
is given via its quadratic form t) £ (u) equal to 

\\duf Xs = J2e d - 1 f (\u'( s ,v)\^^\A Y u\iy ye d S +e d - 2 I |d< dx v 

eeE Jx " Jx " 

where v! is the longitudinal derivative, vl = d s ii 7 and du is the exterior 
derivative of u. Again, t) £ is closed by definition. Adding a potential, we 
define the Hamiltonian H £ as the operator associated with the form t) £ = 
t)(x s ,Q s ) given by 

fie = ||dw-||x e + (u,Q £ u)x c , 

where the potential Q £ is supported in the vertex region X v only. Now we 
use graph result mentioned as an inspiration and choose 

Q £ (x) = -Q(x), 

e 

where Q = Q\ is a fixed bounded and measurable function on X v . The 
reader may wonder that in comparison to (6) the factor e _1 is missing in 
the argument, however, this is due to our choice to perform the squeezing 
by the change of the metric only. 

One can establish the relative (form-)boundedness of H £ with respect 
to the free operator Ax E ■ to a given 77 G (0, 1) there is e v > such that the 
form f) £ is relatively form-bounded with respect to the free form D £ , that is, 
there is C v > such that 

\t) e (u)--0 s {u)\ <riD e {u) + CJu\\ Xe 

whenever < e < e v with explicit constants and . The latter are given 
explicitly in [19]; what is important that they are expressed in terms of the 
parameters of the model which we give below. 

We have mentioned above that our operators acts in different spaces, 
namely the Hilbcrt spaces H — L 2 (G) and H £ = L 2 (X £ ) and their Sobolcv 
counterparts, hence we need to define quasi-unitary operators to relate the 
graph and manifold Hamiltonians. For further purposes we denote 

p e := (voLi_iY;) 1/2 and q(v) = / Qdx v ; 

Jx v 

recall that we introduced the weights p e to be able to treat situations when 
the tube cross sections Y e are mutually different. 

First we define the graph-to-manifold map, J : H — ^ H £ , by 

J/:=e-( d - 1 )/ 2 0(/ e ®T e )®O, (7) 
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where i e is the normalized eigenfunction of Y e associated to the lowest 
(zero) eigenvalue, i.e. ± e (y) = Pe 1 ■ Next introduce the following averaging 
operators 



The map in the opposite direction, J' : H £ — > H, is given by the adjoint, 



In an analogous way one can construct identification maps between the 
Sobolev spaces. They are need in the proofs but not for stating the result, 
hence we refer the reader to [19] for their explicit forms. 

Using these notions in combination with an abstract convergence result 
of [34] one can then arrive at the following conclusions [19]: 

Theorem 4.2. As e — > 0, we have 



for z [Ao,oo). Moreover, <ft(X) = (A — z) _1 can be replaced by any mea- 
surable, bounded function converging to a constant as A — > oo and being 
continuous in a neighborhood ofa(H). 

Corollary 4.1. The spectrum of H £ converges to o~{H) uniformly on any 
finite energy interval, and the same is true for the essential spectra. 

Corollary 4.2. For any A 6 Od; sc (-ff) there exists a family {X £ } with X £ G 
disc(_ff £ ) such that X £ — > A as e — > 7 and moreover, the multiplicity is 
preserved. If X is a simple eigenvalue with normalized eigenfunction <f>, then 
there exists a family of simple normalized eigenf unctions {<ft 6 }s of H £ such 
that \\Jcj) — e ||x E — > holds as e — >• 0. 

4.3. The general case 

So far we have talked for simplicity about the star-shaped graphs only. The 
same technique of "cutting" the graph and the corresponding manifold into 
edge and vertex regions works also in the general case. As a result of the 
analysis performed in Ref. 19 we get 

Theorem 4.3. Assume that G is a metric graph and X £ the corresponding 
approximating manifold. If 




(J'u) e ( S ) = e (d -^ 2 (i e ,u e ( S r))Y. = s 



\\J{H-z)- x -{H e - 
PiH-z)-^ 1 -{H e 



z)-'J\\=0{s 1 /% 
-z)-'\\=0{eV*) 



vol X, 



V 



< oo , sup ||Q \X, 

vev 



inf 4 > 



sup 

vev 



vol dX, t 



v \ I oo 



< OO 



V 



eG-E 
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and 

inf A 2 («) > , inf A 2 (e) > , 

v£V e£E 

where A 2 denotes the second Neumann eigenvalue in the appropriate mani- 
fold region, then the corresponding Hamiltonians H — Aq + q(v)S v and 
H e = Ax E + ^2 v £~ 1 Qv are O(e 1 ^ 2 )-close with the error depending only on 
the above indicated global constants. 

In this way we have managed to solve the problem for quantum graphs with 
r5-couplings under mild uniformity conditions. 

5. Discontinuity at the vertex: the example of 6' s 

While the above results break the Kirchhoff restriction of the previous stud- 
ies, they do not give a full answer; recall that the (5-couplings at a vertex v 
represent a one-parameter subfamily in the n 2 parameter family, n — deg v, 
of all self-adjoint couplings. Let us now investigate the case of S' s as a prime 
example of coupling with functions discontinuous at the vertex. 

5.1. The idea of Cheon and Shigehara 

Our strategy will be the same as before, first we will construct an approx- 
imation on the graph itself and then we we will lift it to the manifold. 
The problem is not easy and its core is the question whether one can ap- 
proximate the ^'-interaction on the line by means of (regular or singular) 
potentials. It was believed for a considerable time that this problem has 
no solution, until Cheon and Shigehara in the seminal paper [9] demon- 
strated a formal approximation by means of of three fS-interaction; a sub- 
sequent mathematical analysis [2,16] showed that it converges in fact in 
norm-resolvent sense. 




Fig. 3. CS approximation scheme on a graph 
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The idea can be extended to ^-coupling on a graph. A scheme of the 
approximation is given on Fig. 3. One starts with a r5-coupling of strength 
6(a) and adds ^-interactions of strength c(a) at the graph edges; the param- 
eter a is the distance of the additional interactions from the vertex. Core of 
the approximation lies in a suitable, a-dependent choice of the interaction 
strengths: we put 

H^ a :=A G + b(a)6 V0 +Y,c(a)5 Vc , b(a) = -A c (a) = ~ 

e 

which corresponds to the quadratic form 
tf' a (f) ■= E ll/ef - ^l/(0)| 2 - \ £ l/e(«)| 2 - domff = W 2 -\G) . 

e e 

Then we have the following result [6]: 

Theorem 5.1. \\{H^ a - z)- 1 - (H fj - z^W = 0{a) holds as a ->■ for 
z ^ R, where || • || is the operator norm on L 2 (G). 

Proof is by a direct computation. We note that the result is highly non- 
generic, both resolvents are strongly singular as a — >• but in the difference 
those singularities cancel. 

5.2. The convergence result 




Fig. 4. Scheme of the lifting 

Now we will lift the above graph approximation result to the manifold 
according to the scheme depicted on Fig. 4. For simplicity assume that 
the star graph in question is finite with all edges having the same length; 
without loss of generality we may put it equal to one. In contrast to the 



December 13, 2010 1:29 



WSPC - Proceedings Trim Size: 9in x 6in delhitalklO 



11 



previous section we have two parameters to deal with, the tube width s and 
the distance of the additional potentials; we choose a — a £ = e a with a € 
(0, 1) to be specified later. The crucial point is the choice of the additional 
potentials. The simplest option is to assume that they are constant, 

e vol X v 

so that J x Q E , v dx — e~ 1 q E (v), where we put 

q s (v ) := b(e a ) = -(3 £ - 2a and q e (v e ) := c(e a ) = -s~ a . 

The corresponding manifold Hamiltonian and the respective quadratic form 
are then given by 

^ = A '.- e " 1 " 28 ^r^- e " 1_B E^- 



eeE 

where \x is the indicator function of the set X, and 

= ||d«||i. -e-^-^IHI^ -e- 1 -^ 

respectively. Note that the unsealed vertex neighborhood X Ve of each of the 
added vertices v e has volume one by construction. 

We employ again the identification operator (7). Using the same tech- 
nique as in the 5 case, one can prove the following result [19]: 

Theorem 5.2. Assume that < a < 1/13, then 

\\(Hi j - l )- l J-J(H^- l y l \\^0 
holds as the radius parameter e — > 0. 

Remark 5.1. The theorem has analogous corollaries as the <5-coupling 
result of the previous section, however, a caveat is due. If < the the 
spectrum of H@' a is uniformly bounded from below as a — > 0. If j3 > 0, on 
the other hand, the spectrum of H^' a is asymptotically unbounded from 
below, inf cr(iJ ,3 ' a ) — > — oo as a — > 0. At the same time, for /3 > the 
spectrum of the approximating operator H@ is asymptotically unbounded 
from below, 'mfa(H^) — > — oo as e — > 0. This fact, existence of eigenvalues 
which escape to — oo in the limit does not contradict the fact that the limit 
operator H 13 is non-negative. Recall that the spectral convergence holds 
only for compact intervals / C K, in particular, cr(H^) n / = implies that 
<j{HP) n I = and aH 13 ' 6 n / = for e > small enough. 

Remark 5.2. While it is easy to see that the parameter a in the approx- 
imation must less than one, the value ^ is certainly not optimal. 
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6. Full solution on the graph level 
6.1. Going beyond S and S' s 

The network approximations of the 8 and (^-couplings described in the 
two previous sections represent the present state of art in this question. 
One naturally asks whether one can extend the technique to other vertex 
couplings. Following the philosophy used here we should look first whether 
such approximations exist on the graph level. 

The simplest extension covers the class of couplings invariant w.r.t. per- 
mutations of edges. It is a two-parameter family containing 5 and 5' s as 
particular cases; in the parametrization (2) its elements are characterized 
by matrices U = a J + bl with \b\ = 1 and \b + adegi;| = 1. The appropri- 
ate approximation in the spirit of Theorem 5.1 was worked out in Ref. 20; 
note that, as with 5 and 5' s , the problem again splits into a one-dimensional 
component in the subspace symmetric over the edges and its complement 
which is trivial from the coupling point of view. 

If we relax the symmetry requirement things become more complicated. 
The first question is what we can achieve by modifications of the original 
Cheon and Shigehara idea, placing a finite number of properly scaled 5- 
interactions on each edge. The answer is given by the following claim [21]: 

Proposition 6.1. Let G be an n-edged star graph and G(d) obtained by 
adding a finite number of 5s at each edge, uniformly in d, at the distances 
0(d) as d — > 0+. Suppose that these approximations yield conditions (1) 
with some A, B as d — > 0. The family which can be obtained in this way 
depends on 2n parameters if n > 2, and on three parameters for n = 2. 

It was demonstrated in Ref. 21 that a family with the maximum number 
of parameters given in the proposition can be indeed constructed. 

In order to get a wider class one has to pass to a more general approx- 
imation. The idea put forward in Ref. 21 was to change locally the graph 
topology by adding new edges in the vicinity of the vertex whose lengths 
shrink to zero in the approximation. This yielded a family of couplings with 
("J ) parameters and real matrices A, B. To get a better result which will 
be described below one has to do two more things: 

• together with adding edges in the vicinity of the vertex one has also 
to remove parts of the graph to optimize locally the approximating 
graph topology, 

• furthermore, one has to add local magnetic fields described by suit- 
able vector potentials to be able to get couplings which arc not 
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invariant w.r.t. time reversion. 
6.2. An alternative unique parametrization 

In order to present the indicated approximation result we have to first 
introduce another form of the boundary conditions (1) derived in Rcf. 7. 

Theorem 6.1. Consider a quantum graph vertex of degree n. If m < n, 
S e C m ' m is a self-adjoint matrix and T £ C m '™~ m ) then the relation 



expresses self-adjoint boundary conditions of the type (1). Conversely, for 
any self- adjoint vertex coupling there is an m < n and a numbering of the 
edges such that the coupling is described by the boundary conditions (8) with 
uniquely given matrices T e £rn,n-m and self-adjoint S G C m,m . 

Remark 6.1. As we have mentioned there are several unique forms of 
the conditions (1). Kuchment [30] splits the boundary value space using 
projections P, Q corresponding to Dirichlet, P^F = 0, and Neumann, Q^S' = 
0, parts and the mixed conditions in the complement. It is easy to see 
that parts singled out correspond to eigenspaces of U corresponding to 
eigenvalues =f1, respectively. The conditions (8) which one call the ST- form 
single out the eigenspace corresponding to — 1 . There is also an analogue of 
(8) symmetric w.r.t. the two singular parts, called PQRS-form, cf. Ref. 8. 

6.3. A general graph approximation 

In view of the above result one can put general self-adjoint boundary con- 
ditions into the form (8) renumbering the edges if necessary. We will now 
describe how those can be approximated by a family of graphs with locally 
changed topology and added magnetic fields. For notational purposes we 
adopt the following convention: the lines of the matrix T are indexed from 
1 to m, the columns are indexed from m + 1 ton. 

The general vertex-coupling approximation, schematically depicted in 
Fig. 5, consists of the following sequence of steps: 

• Take n halflines, each parametrized by x £ R+, with the endpoints 
denoted as Vj, and put a r5-coupling to the edges specified below 
with the parameter Vj(d) at the point Vj for all j = 1, ... , n. 

• Some pairs Vj,Vk, j ^ k, of halfline endpoints are connected by 
edges of length 2d, and the center of each such joining segment is 
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Fig. 5. The scheme of the approximation. All the inner links are of length 2d, some 
may be missing. The grey line symbolizes the vector potential A^ k ^(d). 

denoted as W^j^y. This happens if one of the following conditions 
is satisfied: 

(a) j = 1, . . . , m, k > m + 1, and Tjk 7^ (or j > m + 1, 
k = 1, . . . , m, and Tjy ^ 0), 

(b) j, k = 1, . . . , m, and Sjk 7^ or ( 31 > m + 1 ) 
(T^OAT^O). 

• At each middle-segment point W^^ k } we place a 5 interaction with 
a parameter u^- fc } (<i) . The connecting edges of length 2d are con- 
sidered as consisting of two segments of length d, and on each of 
them the variable runs from zero at W^j^} to d at the points Vj,Vk- 

• On each connecting segment we put a vector potential of constant 
value between the points Vj and 14. We denote its strength between 
the points W^k} and Vj as A^ k ) (d), and between the points W{jj~y 
and Vk as Ar k j\(d). It follows from the continuity that Ar k j\(d) = 
- A (j,k){d) for any pair {j, k}. 

The choice of the dependence of Vj (d), w^^ k } (d) and Aij^\ (d) on the param- 
eter d is crucial for the approximation. In order to describe it we introduce 
the set Nj C {1, . . . , n} containing indices of all the edges that are joined 
to the j-th one by a connecting segment, i.e. 

Nj ={k < m\ S jk ^ 0} U {k < m\ {31 > m + l)(T jt ^ A T u ^ 0)} 

U{k>m + l\T jk ^0} for j < m (9) 
Nj ={k < m\ T kj ^ 0} for j > m + 1 
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We distinguish two cases regarding the indices involved: 

Case I. First we suppose that j = 1, . . . , m and I <E Nj\{l, . . . , m}. Then 
the vector potential strength may be chosen as follows, 

'^argTj, ifRcT, 7 >0, 
^ (arg T jt - tt) if ReT jt <0 
while for vi and to^j} with I > to + 1 we put 

Vl (d) = 1 ~# N l +Jlh=\( T hl) v ^ m + 1) 

= ^ ( — 2+ — | J Vj, Z indicated above, 

where the symbol (•} here has the following meaning: if c e C, then 

f |c| if Rcc > 0, 
{C * ~ \-\c\ if Rec< 0. 

We remark that the choice of Vi{d) is not unique. This is related to the 
fact that for m = rank B < n the number of parameters of the coupling is 
reduced from n 2 to at most n 2 — (n — to) 2 . 

Case II. Suppose next that j = 1, . . . , m and k G iVj D {1, . . . , m}. 

V l=m+l / 

where /x = if 

Re [d-S , J -fc+ ^ >0 
and /x = 1 otherwise. The functions W{j t k} are given by 



w {j,k} = A \2+(d-S jk + J2 TjiTki \ 



and «j (cZ) for j = 1 , . . . , to by 

//AT m / 1 " \ 1 ™ 

« J -(d) = 5 jJ -^2-E( s ifc + d E T * T *)+d E 

fc=l \ l=m+l I l=m+l 

Having constructed the approximating graph we may now investigate 
how the corresponding Hamiltonian behaves in the limit d — > 0. We denote 
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the Hamiltonian of the star graph G with the coupling (8) at the vertex as 
fP tar and jf a P prox w m stand for the approximating operators constructed 
above; the symbols R staT (z) and i?^ pprox (z), respectively, will denote the 
corresponding resolvents. Needless to say, the operators act on different 
spaces: i? star (z) on L 2 (G), while Rf pmx (k 2 ) acts on L 2 (G d ), where G d is 
the Cartesian sum G8(0,d)^"=i^. To compare the resolvents, we identify 
R Ad (z) with the orthogonal sum 

R s d tal (z) = i? star (z)e0, 

which acts as zero on the added edges. Comparing the resolvents is in prin- 
ciple a straightforward task, however, computationally rather demanding. 
Performing it we arrive at the following conclusion [7] which provides us 
with the full answer to our problem on the graph level: 

Theorem 6.2. In the described setting, the operator family /f a PP rox con _ 
verges to if star in the norm-resolvent sense as d — > 0. 

Remark 6.2. There are various modifications of the approximation de- 
scribed above. In Ref. 11, for instance, the S- interactions on the connecting 
segments have been replaced by varying lengths of those segments; the con- 
struction is there performed for scale-invariant vertex couplings, i.e. the 
conditions (8) with S — and any T. 

7. Concluding remarks 

We have demonstrated how one can use scaled Schrodinger operators to 
approximate quantum graph Hamiltonians with different vertex couplings. 
We have worked out the argument for the 6 and (^-couplings. On the graph 
level we have provided a full solution of the problem. 

This suggests how one could proceed further. The approximating graph 
of the previous section has to be replaced by a network with a fat edge 
width e and the r5-couplings by constant potentials of the appropriated 
strength at the segment of fat edge of length e. Similarly the Laplacian is 
to be replaced by magnetic Laplacian on the added edges, the halflength 
of which is set to be d — e a . We call the resulting magnetic Schrodinger 
operator H", where to stands now for the appropriate family of parameters, 
and by H" the corresponding limiting operator on the graph itself. 

Conjecture 7.1. If a > is sufficiently small the approximation result 
analogous to Theorem 5.2 is valid in the described setting for any vertex 
coupling (8) with the same identification operator J. 
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Scaled potentials are not the only way how approximations of nontrivial 
vertex couplings can be constructed. There are other possibilities such as 
replacement of Neumann by suitable position dependent boundary condi- 
tions - for a survey of fresh results we recommend Ref. 35. A more difficult 
question is whether one can accomplish the goal by geometric means. A 
naive inclusion of curvature-induced potentials does not give the answer 
[19] a more elaborate approach has to be sought. 

Let us finally comment on possible physical application of the results 
surveyed here. Thinking of the network as of a model of a semiconductor 
system, one can certainly vary the material parameters. Doping the network 
locally changes the Fermi energy at the spot creating effectively a potential 
well or barrier. From the practical point of view, however, this does not 
help much because our approximations need potentials which get stronger 
with the diminishing tube width e. 

A more promising alternative is to use external fields. In experiment 
with nanosystems one often adds "gates", or local electrodes, to which a 
voltage can be applied. In this way one can produce local potentials fitting 
into our approximation scheme, without material restrictions. This opens 
an rather intriguing possibility of creating quantum graphs with the vertex 
coupling controllable by an experimentalist. 
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